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Abstract. Most differential equations found in chemical reaction networks (CRNs) have the form: 

J = fix) = Svix), 

where a; > 0, that is, x lies in the nonnegative orthant R>o, where 5 is a real d x d' matrix 
{stoichiometric matrix) and u is a column vector consisting of d' real-valued functions having a 
special relationship to 5. Our main interest will be in the Jacobian matrix, f'{x), of f{x), in 
particular in whether or not each entry f'{x)ij has the same sign for all x in the orthant, i.e., the 
Jacobian respects a sign pattern. In other words species Xj always acts on species Xi in an inhibitory 
way or its action is always excitatory. 

In [HKG] we gave necessary and sufficient conditions on the species-reaction graph naturally 
associated to S which guarantee that the Jacobian of the associated CRN has a sign pattern. In 
this paper, given S we give a construction which adds certain rows and columns to S, thereby 
producing a stoichiometric matrix S corresponding to a new CRN with some added species and 
reactions. The Jacobian for this CRN based on S has a sign pattern. The equilibria for the S and the 
S based CRN are in exact one to one correspondence with each equilibrium e for the original CRN 
gotten from an equilibrium e for the new CRN by removing its added species. In our construction 
of a new CRN we are allowed to choose rate constants for the added reactions and if we choose 
them large enough the equilibrium e is locally asymptotically stable if and only if the equilibrium 
e is locally asymptotically stable. Further properties of the construction are shown, such as those 
pertaining to conserved quantities and to how the deficiencies of the two CRNs compare. 



1. Introduction 

In this paper we are concerned with polynomial systems of equations arising from systems of 
ordinary differential equations (ODEs) which act on the nonnegative orthant R^q in W^: 

(1.1) | = /w, 

where / : M^q M'^. The differential equations we address are of a special form found in chemical 
reaction kinetics: 

(1.2) d^ = ^^(^)' 

where 5 is a real d x d' matrix and is a column vector consisting of d' real- valued functions. An 
ODE (jl.ip has reaction form provided it is represented as in (jl.2p with v{x) = [f i • • • Vd'] 
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and 



(1.3) Vj depends exactly on variables Xi for which Sij < 0. 

Since such ODEs are identified with chemical reaction networks, we often refer to these as CRNs. 
Call S the stoichiometric matrix and the entries of v{x) the fluxes. We always assume the 
fluxes are continuously differentiable. Furthermore, in many situations all fluxes Vj{x) are monotone 



dvijx) 



the Jacobian 



nondecreasing in each Xi when the other variables are fixed, that is, v'{x) 

of V, has all entries nonnegative for all x G M^qO. This happens in classical mass action kinetics 
or for Michaelis-Menten-Hill type fluxes. See [BQ, Pa] for an exposition. We shall develop a few 
matrix theoretic phenomena bearing on the properties of f'{x) = Sv'{x), called the Jacobian of 
the reaction. For monotone nondecreasing fluxes the reaction form property (|1.3p is equivalent 
to 



(1.4) 



dvj (x) 



^ 4^ 



dvj{x) 



> and ^ ^ 



Sij < 0. 



dxi \ dxi 

and this is what we shall mostly be using. 

The property analyzed in this paper is whether or not each entry f'{x)ij of the Jacobian of the 
CRN has an unambiguous sign, that is it is the same for all x G K^q. If the sign f'{x)ij is minus 
(resp. plus), then the effect of species j on species i is always inhibitory (resp. excitatory). If such 
is the case we say that /' respects a sign pattern. 

1.1. Sign Pattern of AA^. Given this monotone property, we employ the language of signed 
matrices [BS]. Call a sign pattern a matrix A with entries which are ±aij or 0, where aij are free 
variables. To a real matrix B we can associate its sign pattern A = SP(i?) with itajj or in the 
correct locations. Given a matrix with symbolic entries (i.e., polynomials) we might or might not 
be able to associate a sign pattern. Here, we think of the free variables as being positive. 



Example 1.1. If B 



6 



then A = SF{B) 





-021 



+ai2 

-022 



and 



AA' 



n2 



-012022 Ooi + a: 



012^22 

21 I O22 



Observe that AA^ respects a sign pattern. 



On the other hand, if B 



6 



, then A = SF{B) 



-On 
-021 



+012 
-022 



and 



"^11 + "12 011021 — 012022 
2 2 

011O21 — 012O22 021 + O22 

does not respect a sign pattern. Namely, the off-diagonal entries of AA'' are not positive linear 
combinations of monomials in the Oj. They may attain positive and negative values when evaluated 
at appropriate positive values of the Ojj . 

Theorem 1.2 (cf. [HKG, Theorem 5.1]). Let A be a sign pattern. The hermitian square AA^ of A 
respects a sign pattern if and only if A does not contain a 2 x 2 submatrix whose rows and columns 
can be permuted to obtain a matrix whose sign pattern agrees with the one of 

-1 +1] 



(1.5) 



+1 
-1 



or 



+ 1 +1 
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Such 2x2 matrices either contain 3 minus signs and 1 plus sign, or they contain 1 minus sign and 
3 plus signs. 

Proof of Theorem \1.2[ Suppose the entry {AA'^)ij of AA* fails to respect a sign. As 

(^AA')ij = ^ ^ Aii^Af^j = ^ ^ Aiif^Ajii, 

k k 

this is equivalent to not all terms of the last sum having the same sign. Which is equivalent to 
the existence of k,l with signAjfc = sign74jfc 7^ and sign^j^ = — signAj^ 7^ 0. Hence the 2x2 
submatrix 

■^ik -Ajk 

An Aji_ 

of A will (after a possible permutation of rows and columns) have the same sign pattern as one of 
the matrices in (11.51). ■ 



1.2. Sign pattern for the Jacobian. In the language of chemical reaction networks Theorem 
1.21 has an interpretation as follows: 



Theorem 1.3. The Jacobian f'{x) = Sv'{x) of the right hand side of a reaction form ODE ()1.2p 
with monotone nondecreasing fluxes respects a sign pattern in the positive orthant whenever S does 
not have a 2 x 2 submatrix whose rows and columns can be permuted to obtain a matrix whose sign 
pattern agrees with the one of 

r+1 -1" 
-1 -1 



(1.6) 



Conversely, under mass action kinetics [cf. [BQ,Pa] or Section ^3.2p . if a reversible stoichio- 
metric matrix S contains such a submatrix, then f'{x) = Sv'{x) fails to respect a sign pattern on 
an open dense set of reversible matrices having the same sign pattern as S. 

Since, it is brief we review why this is true. Write S = 5"+ — S- for real matrices 5"+, S- with 
nonnegative entries satisfying the complimentarity property {S+)ij{S-.)ij = 0. If the (z,j)th entry 
of f'{x) = Sv'{x) does not have a sign pattern, then {S^v' {x))ij / and {S^v' {x))ij ^ 0. As 

{S+v'{x))ij = J2iS+)ikv'{x)kj, 

k 

{S-^-v' {x))ij 7^ if and only if for some k, So. > and v'{x)kj ^ 0, so from the reaction form 
property (|1.4p . we get Sjk < 0. To summarize: Sik > and Sjk < 0. Similarly, {S^v'{x))ij 7^ 
if and only if there is some i with {S-)i£ 7^ and v'{x)ij ^ 0, so we get Sj£ < 0. To summarize: 
Sii < and Sji < 0. Taken together this implies that the 2x2 submatrix of S given by rows 
i,j and columns k,i has the same sign pattern as the matrix (jl.6p . up to a permutation of rows 
and columns. The converse reverses this line of reasoning, with the hypothesis requiring robustness 
under small perturbations of S to rule out fluke cancellations. 

See [HKG, §3.1] for details and extensions. 

That many CRNs have Jacobians with sign patterns is the basis for the works of Thomas 
[Th], Kaufman [TK], Soule [So], Gouze [Go], Sontag [AnS, ArS06, ArS08], and many subsequent 
publications. Typically Sontag and collaborators assume this and something considerably stronger 
to obtain results on globally stable equilibria. Sontag has had the philosophy for many years that 
any CRN can be modelled carefully to have Jacobians respecting sign patterns. 
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This paper concerns CRNs whose Jacobians do not have a sign pattern and describes a method 
for transforming such a system of ODEs into a system of ODEs whose Jacobian does have a sign 
pattern, and for which the equihbria of both CRNs remain "the same" (see ^3.21 for a precise 
formulation). We shall refer to this as the sign fixing algorithm. 



2. Fixing the sign pattern for the Jacobian 

In the first part of of this section, §2.11 we describe our sign fixing algorithm. In subsequent 
(sub)sections we show that the algorithm has several (pleasant) properties. Sections ^and ^show 
how equilibria and the classical notion of deficiency behave with respect to our algorithm. 

2.1. An algorithm for eliminating non-signed entries of /'. Let 5 be a stoichiometric matrix 
associated to a chemical network with a submatrix whose sign pattern coincides with that of p.6p . 
which Theorem 11.31 demonstrates is an obstruction for having a sign pattern. Let A, B be the 
species representing the two rows of S corresponding to this bad submatrix. Consider the two 
columns of S belonging to this bad submatrix. 



A 


-Pi 


-P3 


B 


P2 


-Pi 



These yield two reactions in the network of the following form: 

(2.1) piA + Ci^p2B + C2 

(2.2) p^A+p^B + C^^Ci, 

where G N and Cj are some (possibly empty) positive linear combinations of species (avoiding A 
and B). 



We will construct a new network from S to eliminate this bad submatrix. Consider the following 
network, where each reaction of the original CRN remains the same except that we add a new species 
B' , reaction (|2.ip is replaced by 

(2.3) piA + Ci^ B' + C2, 
and we create an additional reaction 

(2.4) B'^p2B. 

Notice that the stoichiometric matrix S associated to this new chemical network will not have the 
bad submatrix we started with. Also, no new bad submatrices have been added in this process. 
Thus we have reduced the number of bad submatrices. 
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S = 





(I23D 


(122]) 















A 











B 





-P4 




P2 













B' 


1 








-1 



We continue applying the same procedure (on this new network) to ehminate any other existing 
bad 2x2 sub matrices. 

In matrix terms, each time we apply this procedure to eliminate a bad submatrix, we change 
one column (e.g. changing reaction (j2.ip into (j2.3p ) of S and we append one additional row (e.g. for 
the "species" B') and column (e.g. for the reaction ()2.4p ) to S. We will call this procedure the sign 
fixing algorithm. 

Definition 2.1. Let S" be a stoichiometric matrix corresponding to a chemical network, and suppose 
S has bad submatrices. We write S for a new stoichiometric matrix with no bad submatrices 
obtained by the sign fixing algorithm applied to each bad submatrix as explained above. S is called 
a sign fixing matrix of 5. 

Each step of the sign fixing algorithm has the interpretation that we are keeping track of ad- 
ditional information. Namely, we measure how much of species B is produced by reaction (j2.ip 
thereby obtaining B' and also we measure how much of species B is consumed by reaction (j2.2p . 
Of course, in the original CRN we just kept track of the net amount of B, a single species, now 
we have replaced this with two species. The added reaction has the eff'ect of identifying the two 
species asymptotically. 

The above definition implies that if S has no bad submatrices, then S = S is the sign fixing 
matrix of itself. If S has multiple bad submatrices, we obtain 5 by a finite number of applications 
of the sign fixing algorithm. 

We emphasize what we have found so far by stating 

Theorem 2.2. Given a GRN one can derive a sign fixed CRN. The Jacobian of the sign fixed GRN 
respects a sign pattern, provided each flux vj is monotone increasing on each Xi . 

This follows from Theorem 11.31 because we have produced a sign fixed CRN S not containing 
the forbidden pattern ()1.6p . 

Example 2.3. We now present a modification of Example [CF05, Table 1.1. (v)] of a CRN who's 
Jacobian fails to have a sign pattern, and we illustrate use of the sign fixing algorithm to create a 
network who's Jacobian has a sign pattern. Consider the network: 

(2.5) A + B ^ F 

(2.6) A + G ^ G 

(2.7) C + D ^ B 

(2.8) G + E ^ 2D 
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The stoichiometric matrix for this network (with species arranged in alphabetical order) is 



S = 



-1 


-1 














-1 





1 


-1 











-1 


-1 


1 


-1 


1 








-1 


1 


2 


-2 














-1 


1 


1 




















1 















The boldface entries are those appearing in bad submatrices of S, of which there are two: 

• Bi, corresponds to species C,D and the forward reactions of (j2.7p and ()2.8p . 

• B2, corresponds to species C,D and the forward reaction of p. 71) and the reverse reaction 
of dSS]). 

The first step of the sign fixing algorithm eliminates B2 and gives 



S = 



The second step eliminates Bi 



S = S = 



-1 


-1 














' 




-1 





1 


-1 
















-1 


-1 


1 


-1 


1 













-1 


1 





-2 


2 
















-1 


1 







1 

























1 































1 





-1 




d gives 














-1 


-1 




















-1 





1 


-1 

















-1 


-1 


1 


-1 








1 








-1 


1 





-2 


2 

















-1 


1 








1 


























1 
































1 





-1 




















1 





-1 



Thus 5 is a sign fixing matrix for S. ■ 

2.2. The story in terms of graphs. To S one often associates a bipartite graph Qs- One set 
of nodes is rows Hi (chemical species) the other set of nodes is columns Cj (reactions). If Sij has 
a — sign (resp. +), then the arrow points from Ri into Cj (resp. out of Cj into Ri). The arrow 
is from Ri into Cj (resp. out of Cj into Ri) if species i is consumed (resp. produced) in reaction 
j, respectively. Which leads us to refer to consumed and produced edges. Also we make a 
distinction between dotted and solid edges, the convention being that all dotted edges touching 
a particular Cj are either all consumed or all produced. Likewise for solid edges. This does not 
determine the choice solid vs. dotted uniquely, for example, completely switching the choice of 
solid and dotted carries the same information. If no reaction is reversible this would be redundant 
information with the direction of arrows, but for reversible reactions dotted vs. solid is needed. 
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This graph is a simphfied version of the species-reaction graph used in [CF06]. Theorem [L2] ( 
also the paragraph fohowing it) in this languages says 



see 



Theorem 2.4. The Jacobian f'{x) = Sv'{x) of the right hand side of a reaction form ODE (jl.2p 
with monotone nondecreasing fluxes respects a sign pattern in the positive orthant whenever the 
graph Qg does not contain a cycle of length four with three consumed edges and one produced edge. 
That is the CRN contains two reactions and two species, one reaction consumes both species while 
the other reaction consumes one species and produces the other. 

The sign fixing algorithm takes a length four cycle C as in the theorem and "breaks it" by 

(1) removing an edge from the cycle 

(2) adding a reaction node C* and species node R* and two edges to the graph 
thereby converting C to "harmless" a cycle of length six. 

Example 2.5. We return to Example 12.31 and observe that its graph is 



B - 



C + D^B 



A + B ^ F 



A 



A + C 



C + E^D 



© 

(Instead of Ri and Cj more descriptive names have been used to denote the nodes.) 

2.3. Uniqueness of the sign fixing matrix. Suppose we have a CRN with stoichiometric matrix 
S. The sign fixing matrix S" of S" is non-unique, since it depends on the indexing of the bad 
submatrices of S. It turns out that this non-uniqueness is easy to classify. One finds that any 
two sign fixing matrices for S can be gotten from the other by "conjugation" with a permutation 
matrix in a certain class which we shall describe completely. 

Suppose B = {bi, . . . ,bn} are the bad submatrices of S. Then S is determined by the finite 
sequence S = Sq, Si, ... , Sn-i, Sn = S, where Sj is the sign fixing matrix of Sj^i with respect to bj. 
Another problem arises if bi and bj share the same positive entry in S. Namely, in this case, Sj will 
equal Sj^i by construction, if i < j. To resolve this problem, we introduce an equivalence relation 
on B: bi ~ bj if and only if bi and bj share a positive entry of S. This is an equivalence relation 
by construction. We will use B to denote the set of all equivalence classes of bad submatrices of 
S, and from now on, we identify each bad submatrix of S with its equivalence class. Thus, the 
number of new columns and rows in S is precisely card(;S). 

Suppose B = {bi, . . . , bn}. Let Sym„ be the symmetric group on {1, 2, . . . , n}, i.e., the set of 
all permutations of n elements. Every a G Sym„ determines a sign fixing matrix S^ of S by the 
finite sequence S, S^^i, . . . , Sa,n-i, Sa,n = S^, where S(j,i is the sign fixing matrix of S with respect 
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to and S^j is the sign fixing matrix of S^j-i with respect to ^^(j)- Clearly, each sign fixing 

matrix of S is determined by a permutation of B, and hence by an element in Sym„. We thus 
identify the set of all sign fixing matrices of S with the elements of Sym^. 

Our result in this subsection gives a map from one sign fixing matrix to another. 

Theorem 2.6. Suppose S is a dxd' stoichiometric matrix corresponding to a CRN with n pairwise 
nonequivalent had suhmatrices. If cr,T E Sym„ then 



'Id 0' 


Sr 


'Id' 


0" 


P 





p 



where P is the n x n permutation matrix associated to t ^cr G Sym^. 

Let C be the matrix obtained by substituting each of the n positive entries of S corresponding 
to the n bad submatrices by 0. By construction, each sign fixing matrix of 5 will have the form 



Sr, 



c 



where No- is anx d' matrix whose rows correspond to the bad submatrices in the order determined 
by a, and Mq- is a d x n matrix whose columns correspond to the bad submatrices in the order 
determined by a. More precisely, the ith row of is the unit vector e'j of length d' if the bad 
submatrix ^^-(j) has a positive entry in column j of S. Similarly, the ith column of is a multiple 
of the unit vector ej of length d if the bad submatrix ^^-(j) has a positive entry in row j of S. The 
multiple is the value of 5 at this positive entry. 

Proof of Theorem 12.61 Note that 

C MrP*' 

P^r -In . 

SO we only need to prove MrP^ = M„ and PNr = No- for the given P. 

Suppose B = . . . , bn}, and e G Sym„ is the identity permutation. With the notation above, 



7 0" 




I 


o" 


t 


P 


Sr 





p 





let 



ai 



Or 



and 



[Pi 



where Oj, i 



I, . . . ,n are the rows of A'g, and /3i, 
a G Sym„, we have 

'""f7(l) 



,n are the columns of Mf. Given 



and Mo- = [Pa 



(1) 



Pa{n)] ■ 



The n X n permutation matrix associated to r ^cr is 

e(r-la)(l) 

P 
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where Cj denotes the unit vector of length n with a one in the ith coordinate and O's elsewhere. By 
construction, PN^- is a matrix whose ith row is the (r~^(T)(i)th row of Nr, so the ith row of PNr 
is a^((^-io.)(j)) = ao-(j). Hence N^j = PN^ as desired. 

To conclude the proof let us verify M^- = M^-P^ . Notice that M^-P^ is a matrix whose ith column 
is the (T~^(T)(i)th column of M^, so the ith column of Mt-P* is = /?cr(i)- This implies 

M„ = MrPK ■ 

3. Equilibria behave well under sign fixing 

Our next goal is to analyze how equilibria for the original CRN compare to equilibria for a sign 
fixed CRN. We shall find that the equilibria are in perfect correspondence. The key to this is a 
simple fact in linear algebra which constitutes the next subsection. 

3.1. Linear algebra associated to sign fixing. We now show how the nullspace of S and the 
nullspace of 5 are related. Likewise for the range of S vs. the range of S. 

Proposition 3.1. Let S be a stoichiometric matrix with a bad submatrix. Let S be obtained from 
S by applying the sign fixing algorithm to eliminate this bad submatrix. Then dim ker S = dim ker 5 
and dim ker S* = dim ker . 

Lndeed, there is a precise correspondence: given v G ker 5 there exists a unique v^q £ K with v = 
[v^ VooY £ ker 5". Conversely, for v = [u* Uoo]* S ker 5, one has v G ker 5. A similar statement 
holds for the left kernels. Furthermore, under this correspondence v has positive {resp. nonnegative) 
entries if and only if v does. 



Proof. Let p and q be the rows of S, and let k and £ be the columns of 5 corresponding to the bad 
submatrix. Without loss of generality assume that the bad submatrix has the same sign pattern 
as [li • For the sake of exposition, here is a picture: 



k I 




S 







k e 




d' + 1 













p 











q 






















d+l 





1 





-1 



Suppose S E M'^^'^' and let v G ker 5. Then clearly v = [f* ve]^ £ kerS. Moreover, if t; G 1K>0) 

then V € ^"^q^- Conversely, every v = [f* Voo]* S ker 5 satisfies ve = v^o and so gives rise to 
V S ker S. Again, positivity is preserved. The corresponding analogous statements and proofs for 
the left kernel are left as an exercise for the reader. ■ 

Remark 3.2. After applying one step of the sign fixing algorithm the number of bad submatrices 
decreases. More precisely, let 5" be a stoichiometric matrix with a bad submatrix. Let S be obtained 
from S by applying the sign fixing algorithm to eliminate this bad submatrix. Then the number of 
bad submatrices in S is less than the number of those in S. 
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Indeed, the form of S, because the added row and column have all entries but two equal to zero, 
guarantees S does not contain the one bad submatrix under attack, and at the same time no new 
bad submatrices have been added. Thus we have reduced the number of bad submatrices by at 
least one. ■ 

Remark 3.3. The sign fixing algorithm for the situation of the hermitian square AA^ of a sign 
pattern A, goes just as in ^2.11 with A replacing S in the picture in the proof of Proposition 13. 1[ 

Details are left as an exercise for the interested reader. ■ 



3.2. Behavior of equilibria and steady states in mass action kinetics under sign fixing. 

This subsection uses the linear algebra result of the previous subsection to show that the equilibria 
of a CRN and of its sign fixed CRN are in perfect correspondence. We shall show this for mass 
action kinetics, although as one will see from the arguments here it works for a much more general 
class of CRNs. 

We now review mass action kinetics with the primary aim of introducing our notation. The 
postulate of mass action kinetics is "the reaction rate is proportional to reactant concentrations" . 
For instance, for the chemical reaction 

2A + B ^AC 

the reaction rate is k2A+B^4cx\xB, where x denotes the concentration of a species and A;2a+b^4c > 
is the rate constant. The corresponding ODE is 



XA 




XB 




Xc_ 





In general, for S £ 



— 2k2A+ 3-^40 ^A'^B 

the flux vector v{x) is given by 





'-2 




1 




4 



x. 



' min{0, Sji} 



J. 



(Here /cj > is the rate constant associated to the ith reaction and Xj is the concentration of the 
jth species.) 

If the ODE (jl.ip admits a positive vector in the left kernel (i.e., there exists m G M^g with 
m ■ X = m ■ f{x) = 0), then the ODE is called conserving. This reflects quantities (like the mass 
or the number of carbon atoms) being conserved. An obvious sufficient condition for ODEs of the 
form (jl.2p is ker 5"* H M^g ^ {0}. If this is satisfied, we say that S is conserving. By Proposition 
13.11 this condition is preserved under the sign fixing algorithm. 

Corollary 3.4. Let S he a stoichiometric matrix, and suppose S has a had submatrix. Let S 
he obtained from S hy applying the sign fixing algorithm to eliminate this bad submatrix. If S is 
conserving, then so is S. Moreover, the reaction form differential equations (jl.2p associated under 
mass action kinetics to S and to S, respectively, have the same equilibria in the following sense. 



Suppose S G M''^'* . // X 
Sv{x) = 0. Conversely, if x £ 



I'^Q {resp. X G I 



(resp. X G M^"]^^) satisfies Sv{x) = 0, then 
i>g) satisfies Sv{x) = 0, then there exists a unique 
joY . {Here v will be used to denote a 



Xoo £ ]R>o {resp. Xoo S M>o) with Sv{x) = for x = [x* 
flux vector associated under mass action kinetics to S.) 
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Proof. This is essentially a consequence of Proposition 13.11 and the figure contained in its proof 
describes the notation we now use. Let p and q be the rows of S, and let k and i be the columns 
of S corresponding to the bad submatrix, and assume without loss of generality Sq£ > 0. 

If v{x) G kerS'nM^^^ and x = [x* Xoo]*, then (by construction) the first d' entries of v{x) 
coincide with v{x), that is, 

v{x)i = v{x)i, i = l,...,d'. 

Additionally, 

(3.1) = Xd+i = v{x)i - v{x)d'+i 
so we obtain 

(3.2) v{x)d'+i = v{x)i = v{x)i. 

Note that by construction, v{x)ii'+i depends only on x^o and thus we can solve (j3.2p for x^o uniquely. 
Hence 

d'+l d' 

= ^9 = X] ^qi^{^)i = ^^SqiV{x)i + Sqiv{x)d' + l 
i=l i=l 



SqiV{x)i + Sqiv{x)i = ^ SqiV{x)i. 

i = i 1=1 



For s ^ q, 



d'+l d' 

= = ^ Ssiv{x)i = y^^Ssiv{x)i 

i=l i=l 

proving Sv{x) = 0. (Alternatively, the conclusion can be reached by the proof of Proposition 13.11 1 
The calculation above reverses to show that converses of these implications hold as well. ■ 

Theorem 3.5. Let S he a stoichiometric matrix corresponding to a chemical reaction network. 
Then the reaction form differential equations corresponding to S and to its sign fixing matrix S 
under mass action kinetics have the equilibria which are equivalent under the correspondence in 
Corollary 13.41 



Proof. This follows easily from Corollary 13.41 and Remark 13.21 by an induction on the number of 
bad submatrices of S*. ■ 

Remark 3.6. Theorem 13. 51 and Corollary 13.41 extend to more general, reaction form ODEs (jl.2p with 
monotone fluxes. In one step of the algorithm the key is to add a reaction consuming exactly one 
(new) species (variable) x^o- Since this is an artificial reaction we can specify a flux v{x)d'+i and 
the key is to pick it to be monotone and surjective, e.g. it depends only on Xoo and is linear. This 
ensures the solvability of ()3.2p for Xqo • The uniqueness of Xoo is then guaranteed by the monotone 
property. Under these assumptions both proofs work verbatim. ■ 



12 



HELTON, KATSNELSON, AND KLEP 



3.3. Local stability is preserved by sign fixing. In the previous subsection we showed that 
the equilibria of the original CRN sit in a perfect correspondence with those of the sign fixed 
CRN. An important question is whether or not stability of an equilibrium of the original CRN 
implies stability of the corresponding equilibrium of the sign fixed CRN. This question is open to 
interpretation because the sign fixing CRN contains a rate constant which we are allowed to define. 
Let us call this rate constant k. A natural version of the question would be: is there an a priori 
choice of k such that the equilibrium of the original CRN is stable if and only if the corresponding 
equilibrium is stable for the sign fixed CRN. While we have not analyzed global stability, we have 
analyzed and answered the question for local asymptotic stability. We found that if we choose k 
large enough, then one of the eigenvalues of the sign fixed Jacobian will be very negative, and all 
the others will be close to the eigenvalues of the Jacobian of the original CRN. Recall that a matrix 
is said to be stable if all its eigenvalues have negative real part. An equilibrium xq G M>q of an 
ODE of the form (jl.ip is locally asymptotically stable if the matrix /'{xq) is stable. 

As before, we assume mass action kinetics although this assumption can be weakened to reaction 
form ODEs ()1.2p with monotone entrywise surjective fiuxes. 

Theorem 3.7. Let S be a d x d' stoichiometric matrix with a bad submatrix. Let S be obtained 
from S by applying the sign fixing algorithm to eliminate this bad submatrix. Write J{x) = Sv'{x) 
and Jfc(x) = Sv'{x). Here k denotes the rate constant assigned to the additional reaction created 
in the sign fixing algorithm. Fix a point x G K>^^ and let x G M>q denote its first d components. 
Furthermore, let J = J{x) and Jk = Jk{x). 

Then d of the eigenvalues of the {d + 1) x (d + 1) Jacobian matrix Jk (counting multiplicity) 
converge (as k — > oo) to the d eigenvalues of J and the remaining eigenvalue is real and converges 
to — oo. 



Without loss of generality, we may assume the bad submatrix in S is the 2x2 bottom right block 



and Sd^d' > 0. Then the relationship between the dx d matrix J 
matrix Jk = Sv'{x) is as follows: 



Sv'{x) and the {d+l)x {d+l) 



Jk 



Jd~l,l 

Jd,i - Sd,d'-5^ 



J 







dxi 









Jd- 



l,d 













+ 


— Sd4' 




1 



dxd 



9'"di (^) 
dxi 



Jd,d Sd4' -Q^ 





kSd4' 
-k 



d^d' (^) 
dxd 



Let c(A) = det(J — Xld) G M[A] and Cfc(A) = det(Jfc — Xld+i) G ]R[A] denote the characteristic 
polynomials (in A) of J and Jk, respectively. 

Lemma 3.8. The degree d+l polynomials ^c^ converge uniformly on compact subsets of C to the 
degree d polynomial — c. 
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Proof. Let us consider J^ — XI and its determinant. For notational convenience let us write s = S^^' 
and Vd'j = Then 



Ck = det( Jfc - M) 





Ji,i - A 




J 


U 




= det 


'-'a— 1,1 




h 1 J 


n 






Jd,l - SVd',! ■ 


Jd,d ~ 


- SVd' ,d — 


X ks 






Vd',1 




Vd',d 


-k-X 






"Ji,i - A • • • 


Jl,d 









= det 


Jd-l,l 


Jd~l,d 











Jd,l 


Jd,d — A 


-sX 








Vd',1 


Vd',d 


-k-X 










'Ji,i - A 




Jl,d-l 


Jl,d 


= {-k 


— A)c + sX det 










Jd-l,l 


■ ■ ■ Jd-l,d-l - A 


Jd-l,d 






Vd',1 




Vd',d-1 


Vd',d _ 















Thus 

(3.3) det( Jfc - XI) = Ck = {-k - X)c + sXh = {-k - A) det J + sXh. 

Note h G K[A] is a polynomial of degree < d — 1 in A and does not contain k. Thus 



1 



-A; - A 



-Ck 



k k 
uniformly on compact subsets of C. 



1 . , fc— »oo 

—sXn > 

k 



In fact, the polynomial h from the proof of Lemma [3. 81 is of degree < d — 2. Since the ODEs are 
in reaction form, Sd,d' > implies Vd',d = §^ = 0, cf. (II. 4p . 

Proof of Theorem 13.71 Let Xj , nij denote the zeroes of Cfc together with their multiplicities and 
Xj, rrij denote the zeroes of c. Certainly Cfc is analytic in the complex variable A, thus d zeroes of Cfc 
(counting multiplicity) converge to the zeroes of c. This is a standard consequence of the argument 
principle, since we can put a small circle Cs around a zero of Xj and for large enough k the winding 
number (with respect to 0) of ^Cfc on Cg equals that of c. Thus c and Cfc have the same number of 
zeroes inside Cg. 

Similarly, to analyze the point at infinity, one can draw a circle Cpt of arbitrarily large radius R 
containing all zeroes of c. The winding number (with respect to 0) of c around R is d, so for large 
enough k the winding number of ^Cfc is also d, thus one zero of Cfc, without loss of generality denote 
it ^d+i ^^^^ outside of Cr. Hence the sequence x^^^ diverges to infinity. Since all coefficients of the 
polynomial Ck are real, its zeroes are either real or occur in conjugate pairs. So x^j^^ must be real, 
since if not Cfc would have two zeroes outside of Cr. 

Let us retain the notation from the proof of Lemma 13. 8[ Then 

Cfc = {-k- X)c + sXh e R[A], 
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where c G ]R[A] is of degree d, h & IR[A] is of degree < d — 1, and s G M. Thus for A > big enough, 
—Ac dominates sXh. For such A > the sign of Cfc(A) wih equal the sign of — c(A) for any A; > 0. 
This shows that with R big, the zero x^_^^ of outside of Cr must be negative, thereby concluding 
the proof. ■ 

Theorem 3.9. Let S be a stoichiometric matrix corresponding to a CRN and let S be one of its 
sign fixing matrices. Then there exists a choice of rate constants for the added reactions such that 
the equilibria of the reaction form ODEs corresponding to S are locally asymptotically stable if and 
only if the same holds for the equilibria of the reaction form ODEs corresponding to S. 

Proof. This follows easily from Theorem 13.71 by an induction on the number of bad submatrices in 
S. m 

3.4. Craciun-Feinberg theory: determinants of Jacobians. This brief subsection is for those 
familiar with the Craciun-Feinberg (CF) theory [CF05, CF06] and we observe that it behaves well 
under sign fixing. Recall the key requirement of the CF theory in order to invoke its consequences is 
that the determinant A of the Jacobian has an unambiguous sign. Note from (j3.3p with A = that 
the determinant of the Jacobian of the sign fixed CRN is a scalar times A. Thus one determinant 
has an unambiguous sign if and only if the other one does. 

The conclusion is that if the CF theory applies to a CRN, then it applies to the closely related 
CRN whose Jacobian respects a sign pattern. 



4. Deficiency vs. sign patterns 



An important notion in chemical networks is that of deficiency. In this section we show that 
sign fixing might increase the deficiency of a CRN by at the most the number of bad submatrices 
for the stoichiometric matrix of the original CRN. 

We follow the notation and terminology of Gunawardena [Gu] (or see [Fe, CDSS]). Thus, we 
denote: 

n := the number of complexes of the network, 
i := the number of linkage classes of the network, 
s := the rank of the stoichiometric matrix, 
and the topological deficiency of the network is 

6 := n — £ — s. 

4.1. Zero deficiency vs. sign patterns. A natural question is whether the sign pattern of f'{x) = 
Sv'{x) has any correlation to the CRN having zero deficiency. The answer is no, and in this 
subsection we give examples of 

(1) chemical networks S with zero deficiency and no sign pattern for Sv'{x); 

(2) chemical networks S with nonzero deficiency and a sign pattern for Sv'{x). 

Example 4.1 (see [Ka, §4.3] for more details). Consider the reaction network 

B 

B^C 
C^A + B. 
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The deficiency of the CRN is easily seen to be zero. However, 



S 



-1 



1 




1 



1 



1 
1 

-1 



so Sv'{x) will not respect a sign pattern (Theorem II. 2p . It will have exactly one entry without 
a sign. In order to obtain an example of a deficiency zero network with an arbitrary number of 
non-signed entries in the Jacobian, one simply considers a network with the following stoichiometry: 

^5 



Conversely, having a sign pattern will not yield any information about the deficiency of the 
network. 

Example 4.2. Consider the chemical reaction network 



B^B' ^C^C 



with stoichiometric matrix 



-1 


-1 











1 


1 











1 





-1 








-1 





1 











1 


1 


-1 








-1 


-1 


1 





1 








1 


-1 


-1 








-1 


1 





1 








1 





-1 








-1 



By Theorem 11.21 Sv'{x) respects an unambiguous sign pattern but the deficiency of the network is 
one. To achieve arbitrary deficiency one can employ a block diagonal construction as above. 

4.2. Deficiency and the sign fixing algorithm. In this subsection we consider how the defi- 
ciency of a CRN changes after we apply the sign fixing algorithm to produce a new CRN. 

Let Si be the stoichiometric matrix for a CRN, and let 5*2 be the sign fixing matrix of 
with respect to some bad submatrix. All variables with subscript 1 refer to the original CRN and 
variables with subscript 2 refer to the new CRN unless otherwise noted. Also, C denotes the set 
of all complexes of a network and C denotes the set of all linkage classes of a network. We also 
denote A5 := 82 — 81, An := ?i2 — rii, Ai := £2 — ii, and As = S2 — si. Assuming that Si has a 
bad submatrix corresponding to the species A, B, then this network has 2 reactions of the form 



(4.1) 
(4.2) 



PiA + Ci 
P3A + piB + C3 



P2B + C2 
C4, 



where pi,p2,P3,P4 G N and Ci,C2,C3,C4 are some (possibly empty) positive linear combination 
of species. The only changes to the new network are we add a new species B', reaction ()4.ip is 
replaced by 



(4.3) 



PiA + Ci^B' + C2 
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and we create an additional reaction 
(4.4) 



B' ^P2B. 



By Proposition 13. H we always have As = 1. To get a better handle on the change of deficiency, we 
proceed as follows. 

Lemma 4.3. If Si is the stoichiometric matrix of a CRN and S2 is its sign fixing matrix with 
respect to a bad submatrix, then the following inequalities are sharp: 



(4.5) 



A£<2 andl< An< 3. 



Proof. With the notation above, the only possible new complexes are B' + C2, -B', and P2B. Hence 
An < 3. Also, the only possibly new linkage classes are [B' + C2]2 and [p2-B]2, so Ai < 2. The fact 
that An > 1 is obvious. 



To show that the inequalities (|4.5p are sharp, consider the following network: 

B + 2C ^5D 
A + B^C. 

This CRN has the stoichiometric matrix 



Si 



Hence ni = 5 and ii = 2. The sign fixing matrix for 5i (with respect to the species A, B) is 



-1 


-1 





1 


-1 


-1 


2 


1 


-2 








5 



-1 


-1 











-1 


-1 


1 


2 


1 


-2 











5 





1 








-1 



S2 



and our new chemical network is 

A^2C + B', A + B 
B + 2C ^ 5D, B' B. 

We thus see that n2 = 8 and £2 = 4, so An = 3 and Ai = 2. ■ 

We now need an efficient way to determine which complexes and reactions of the new CRN 
affect An and Ai. We shall define functions which will precisely determine which complexes and 
linkage classes in the new CRN increase An and Ai. Define (j) : {B' + C2, B',p2B} {0, 1} by 



(t>{B' + C2) 
<j>{p2B) 



2 if C2 / and P2B + C72 G C2 

1 otherwise 

'1 iip2B^Ci 

otherwise 
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We also define ip : {[B' + C2]2, [B'h} ^ {0, 1} such that 



1 if P2B + C2£C2 and [B' + C2]2 R [P2B + C2]2 
otherwise 



1 if P2B Ci and C2 / 
otherwise 



The advantage of this new notation is that we now have a succinct way to measure An and A£ : 

(4.6) An = 4>{B' + C2) + <I){P2B) 

(4.7) Al = i^{[B' + C2]2)+^{[B']2). 

Equation (j4.7|) follows directly from the definition of AH. and the construction of the new network. 
However, ()4.6p needs more justification. 

Lemma 4.4. VFii/i the setup described above, (j4.6p holds. 



Proof. By construction, notice that we always have Ci C C2 or Ci \ {P2B + C2} C C2. Also, 
{B' +C2, B' ,p2B} are the only possible new complexes that are not inCi. Observe that ifp2B+C2 S 
C2, then C2 = Ci \J{B' + C2, B',p2B}, whence 

An = card(Ci U{S' + C2,5',P25}) -card(Ci) 

= (card(Ci U {B' + C2, B'}) - card(Ci)) + (card(Ci U {P2B}) - card(Ci)) 
= card(Ci U {B' + C2, B'}) - card(Ci) + <^(p2S), 
where the last equality follows by the definition of (p. Also, if C2 = 0, then 

card(Ci U {B' + C2, B'}) - card(Ci) = 1; otherwise, card(Ci U {B' + C2, B'}) - card(Ci) = 2, 

so (j){B' + C2) = card(CiU{S' + C2, S'})-card(Ci) by construction. This implies An = (t>{B' + C2) + 
(j){p2B), as desired. On the other hand, suppose p2^+ (^2 C2. If C2 = 0, thenp2-B+C2 = P2B G C2, 
which is a contradiction. Thus we must have C2 7^ 0. A simple count shows 



2 if P2B Ci 
1 otherwise 



An = card(Ci \ {P2B + C2} U {B' + C2, B',p2B}) - card(Ci) 

= cPiB' + C2) + HP2B), 
where the last equality follows directly from the definition of (j). ■ 

Theorem 4.5. Let Si be the stoichiometric matrix to a chemical network with a bad submatrix, 
and let S2 be the sign fixing matrix with respect to this bad submatrix. Then < A5 < 1, and this 
inequality is sharp. 

Proof. First, Lemma shows A(5 = An - A^ - 1 < 3 - A£ - 1 < 2. Notice that if A6 = 2, 
then An = 3 and A£ = 0, so cl){B' + C2) = 2, (^(^2-6) = 1 and V([5' + C2]2) = V'([-B']2) = 0. 
However, (j){B' + C2) = 2 implies C2 ^ 0, and (j){p2B) = 1 implies P2B Ci by construction. Hence 
ip{[B']2) = 1 by the definition of tjj and this is a contradiction. Thus, we cannot have An = 3 and 
A^ = 0. This proves A6 < I. 

Now suppose A5 < to derive a contradiction. By Lemma 14.31 An = 2 and Ai = 2, or An = 1 
and A^ > 1. 
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Case 1: Suppose An = 2 and A£ = 2. Notice if A£ = 2 then ^([S' + C2]2) = i^UB'h) = 1 and 
thus C2 7^ 0, P2B + C2 G C2, andp2-B Ci by construction. Hence 4){p2B) = 1 and 4>{p2B + C2) = 2 
by construction, so An = 2 + 1 = 3 by (|4.6p and this is a contradiction. 

Case 2: Suppose An = 1 and A^ > 1. Notice, An = 1 and (jiTBl) imply (/>(p2^) = and 
P2^ G Ci, so il){[B']2) = and + 02)2) = 1 since M > I implies p2-B + C2 G C2 and 

[B' + C2]2 n [p2B + C2]2 = 0. If C2 = 0, then [B' + C2]2 = [5']2 = [^2^12 = [P2B + C2]2. Hence 
[B' + C2]2 n \p2B + C2]2 / and this is a contradiction. Thus C2 + 0. But then (\>{B' + C2) = 2 
by construction, and An > 2 by (j4.6|) : contradiction. 

To show that these inequalities are sharp, consider the following network: 

2A ^ 3B + C 
A + B C. 

This CRN has the stoichiometric matrix 

" -2 -1 2 

Si= 3 -1 - 

1 1 - 



Hence ni = 4 and ii = 2. The sign fixing matrix for Si with respect to A, B is 



S2 



-2 


1 


2 








-1 


-3 


3 


1 


1 


-1 





1 








-1 



and our new chemical network is 



We thus see that n2 = 7 and £2 
3-1-1 = 1. 



3B + C 2A^B' + C 
A + B C 
B' 3B. 

= 3. So An = 3 and Ai = 1. Thus, A6 = An - Ai - As 



Corollary 4.6. Suppose S is the stoichiometric matrix for some chemical network, and it contains 
k bad submatrices. If S is its sign fixing matrix, then < A6 < k. 

Proof. Recall that by definition, S is determined by a recursive sequence of at most k sign fixing 
matrices, each with respect to a certain bad submatrix from the previous matrix in the sequence. 
An application of Theorem 14. 5 1 at each step yields our desired result. ■ 



In fact, the upper bound for A5 in Corollary 14.61 is the number of equivalence classes of bad 
submatrices of S as defined in ^2.31 Also, all sign fixed matrices S obtained from our algorithm 
have the same deficiency (independent of the order in which the sign fixing algorithm is applied). 

Theorem 14.51 together with formulas ()4.6p and (j4.7|) helps determining necessary conditions 
for A6 = 1 for a stoichiometric matrix and its sign fixing matrix with respect to a certain bad 
submatrix. Notice that A6 = 1 implies An = 2 and AI = or An = 3 and AI =1. If C2 = 0, 
then P2B = P2B + C2 G Ci, so (t){p2B) = 0. Also, C2 = implies (j){B' + C2) = 1 by construction. 
So An = 1 by (14. 6p . and hence, A6 = by Theorem 14.51 This observation yields the following: 



SIGN PATTERNS FOR CHEMICAL REACTION NETWORKS 



19 



Theorem 4.7. Let S be the stoichiometric matrix for a chemical network. Suppose that the column 
corresponding to each bad submatrix of S with the positive entry has only one positive entry. Then 
if S is a sign fixing matrix for S, we have AS = 0. 

Proof. By assumption, each bad submatrix of S corresponds to 2 reactions of the form: 

(4.8) piA + Ci ^ p2B 

(4.9) p3A+piB + C2^C3, 

where A, B are species, Pi,P2,P3,P4: G N and ci,C2,C3 are some (possibly empty) positive hnear 
combination of species. As aheady shown, the deficiency for the sign fixing matrix of S with respect 
to this bad submatrix does not change. An inductive procedure yields our desired result. ■ 



5. An ALTERNATIVE SIGN FIXING ALGORITHM? 



-2 


-1 


4 


4 


-4 


-12 


4 


4 








4 


-1 


-2 








10 


-2 


-6 


-4 


4 



Given a stoichiometric matrix S with bad submatrices, there is an easier way of eliminating 
these. Instead of performing the sign fixing algorithm for each submatrix separately and thus 
adding a row and a column in every step, we can add only one row and column and eliminate all 
bad submatrices in a single step. Unfortunately, this construction changes the dimension of ker S; 
thus the two matrices yield reaction networks with very different equilibria structure. We illustrate 
this with an example. 

Example 5.1. Suppose 



S 



Notice that S has several bad submatrices. We start by adding a row and column of zeros to S. 
Pick a bad 2x2 submatrix of S. Replace the positive entry Spi of S by 0, add +1 to the ^th 
entry of the new row and add Spi to the pth entry of the new column. Repeat this for all the 
bad submatrices. After all the bad submatrices have been eliminated, the bottom right entry is 
changed into the negative sum of all the entries in the last row. We obtain a matrix S with no bad 
submatrices. In our example this is 



S 



We note that 



-2 


-1 








-4 


8 


-12 





4 








4 





-1 


-2 








4 





-2 


-6 


-4 





14 


1 


1 


1 


1 


1 


-5 







1 


1]'. 


[1 2 


1 



and 



ker 5* = span I [1 1 1 l]*}. 
The corresponding kernels for S are: 

ker 5 = span I [2 4 1 3 2]* 



} 
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and 

ker5* = {0}. 

Given its kernel, the ODE x = Sv{x) can have equihbria only on the boundary. In fact, each 
solution to Sv{x) = can be shown to satisfy v{x) = 0. 

Also, assuming mass action kinetics, 

v{x) = [kixfxl"^ k2XiX3x\ k^x'^xf k^xf k^xf]* 

so for every X2 G ^>o, 

u3/2. /T- .4/7-^18 



_ k2\^ _ Ak'l'^^/k^^xl^ _ k2 

■^1 7t= n — An — ' -^3 , o An — ' -^4 



2^/klVh\^x^' kl^ ' 2^1^x1 

yields a positive solution to Sv{x) = 0. 

Hence it is not possible to recover positive equilibria for the chemical CRN by S from those 
obtained by S. 

Also note that there is a nonnegative vector orthogonal to the range of S*, thus the corresponding 
reaction form dynamics has a conserved quantity. On the other hand, S is not conserving. 



6. Software 

The discovering of the results in this paper was considerably facilitated by computer experiments. 
The programs we wrote to do this might be of value to a broad community, so we documented 
them and provided tutorial examples. They are found on the web site 

http : //www . math . ucsd . edu/~cheiiicomp/ 

The Mathematica files provided contain software for dealing with equations that come from 
CRNs; dx/dt = f{x) = Sv{x) as in (jl.2p . Some of our commands focus on the Jacobian, /', of /; 
they do the following 

(1) compute the Jacobian /' of / (given say the stoichiometric matrix S); 

(2) check existence of a sign pattern for f'{x) which remains unchanged for all x > 0, using Theorem 
11.21 in this paper; 

(3) implement the sign fixing algorithm in §2.1t 

(4) compute the Craciun-Feinberg (CF) determinant [CF05, CF06, CTF] of /' (governs CRNs with 
outflows for all species with outflow rate constants equal to one); 

(5) compute the more general Helton-Klep-Gomez core determinant [HKG] of /' (governs CRNs 
with any number of outflows). 

The CF determinant and core determinants are used in tests to count the number of positive 
equilibria, namely x* > such that f{x*) = 0. (For more information, please look at the papers 
[CHW] and [HKG] and the original Craciun, Feinberg et al. papers, [CF05, CF06, CTF].) 

Another part of our Mathematica package deals with deficiency of reaction form differential 
equations, as discussed in 21 The software allows us to compute the deficiency of a CRN as 
well as conversion of representations as follows. One starts with the traditional representation 
f{x) = Sv{x). Our program produces the representation 



Sv{x) = YAk^p{x) 
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where is the Laplacian of the "complexes graph" of the chemical reaction network. Y is the 
matrix whose columns are indexed by complexes and which contain nonzero entries corresponding 
to chemical species which enter the complex. ■0 is a list of monomials in the chemical concentrations. 
For details, see [Gu,Fe,HJ]. 

Our commands also compute the components of the complexes graph. Capability to automati- 
cally plot planar graphs is under development and should be available soon. 

7. Conclusions 

(1) The effect of any species on another species is always (for all concentrations) inhibitory or 
always excitatory, if (generically only if for reversible CRNs) the species-reaction graph of S 
has no bad cycle. 

(2) If a CRN has bad cycles, then the sign fixing algorithm produces a CRN whose Jacobian respects 
a sign pattern. The key properties are: 

(a) Equilibria in perfect correspondence; 

(b) Local asymptotic stability corresponds perfectly; 

(c) Conservation laws correspond; 

(d) The determinant of the sign fixed Jacobian is a constant multiple of the original Jacobian, 
thus the theory of Craciun and Feinberg applies to both CRNs or to neither; 

(e) Deficiency does not drop and its increase has a simple bound. 

Possibly the tests and constructions outlined here will be useful to some experimentalists. In 
an experiment where one aims to understand which reactions occur and the pattern of inhibition 
and excitation, one "first" obtains the species-reaction graph Q. Based on measuring some limited 
number of concentrations one determines the inhibitory or excitatory effects. An issue is whether 
or not enough concentrations were measured, for a possibility is that species Xi has an excitatory 
effect at some concentrations and an inhibitory effect at others. 

The results here give a simple way to sort out this problem. First suppose the CRN does not 
contain two reactions involving two species A and B, one species being consumed by both reactions 
while the other species is consumed by one reaction and produced by the other. Then there is 
no such problem (see Theorem \1.3^ : in principle, one concentration measurement (of all species) 
suffices. If G has such a bad pair of reactions, probably there is trouble. Our sign fixing algorithm 
shows additional measurements which, if they can be made, fix this trouble. 

We emphasize that our analysis does not say what effect a species i has on species j in a single 
reaction but it bears on its effect inside the entirety of the CRN. 
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